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ABSTRACT 

We present the first observations of the circumstellar (CS) disk system 51 Ophiuchi with the Far 
Ultraviolet Spectroscopic Explorer. We detect several absorption lines arising from the unusual metastablc 
atomic species N I ( 2 D), N I ( 2 P), and S II ( 2 D). These levels lie 1.8 - 3.6 eV above the ground level and 
have radiative decay lifetimes of 2 days or less, indicating that the lines arise from warm CS gas. The high 
S/N FUSE spectra, obtained six days apart, also show time- variable absorption features arising from N I, 
N II, O I and Fe III, which are redshiftcd with respect to the stellar velocity. The resolved redshifted 

absorption extends over many tens of km s _1 (40 for N I, 100 for N II, 65 for O I and 84 for Fe III). 

We calculate column densities for all the variable infalling CS gasses, using the apparent optical depth 
method. The Fe III and N II infalling gasses must be produced through collisional ionization, and the 
ionization fraction of nitrogen suggests a gas temperature between 20000 and 34000 K. The infalling gas 
shows a peculiar, non-solar composition, with nitrogen and iron more abundant than carbon. We also 
set upper limits on the line-of-sight column densities of H2 and CO. These observations strengthen the 
connection between 51 Oph and the older debris-disk system Pictoris, and indicate that there may be 
infalling planctesimals in the 51 Oph system. 

Subject headings: stars: individual (51 Ophiuchi) — circumstellar matter — planetary systems: 
formation — comets: general 
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1. INTRODUCTION 

51 Ophiuchi (HD 158643) is a young B9.5 Ve star 
(~3x 10 5 yrs old; van den Ancker, de Winter, & Tjin A 
Djie (1998)), showing a large infrared emission excess aris- 
ing from circumstellar (CS) dust. This dust is presumed to 
be located in a disk similar to that observed around (3 Pic- 
toris (Waters et al. 1988). However, although 51 Oph lies 
only 131^3 pc away from the Sun (van den Ancker, de 
Winter, & Tjin A Djie 1998), its CS dust disk has not 
yet been imaged. ISO spectroscopy of the 51 Oph CS 
dust shows the presence of crystalline silicates, which are 
seen in the spectra of (3 Pic and Solar System meteorites 
and comets, but not in spectra of interstellar (IS) dust 
(Waelkens et al. 1996). The disk of 51 Oph is probably 
close to edge-on, because the very large projected rota- 
tional velocity of the star (v smi — 270 km s _1 ; Dunkin, 
Barlow, & Ryan (1997)) suggests that our line-of-sight is 
parallel to the equatorial plane of the system. 

51 Oph also has a CS gas disk, as evidenced by double- 
peaked hydrogen emission (Slettebak 1982) and narrow, 
non-IS absorption features in the ultra-high-resolution 
spectrum of the star (Crawford, Craig, & Welsh 1997). 
Cold CS C I gas (T rs 20 K) has also been observed to- 
ward 51 Oph; this species has a very short photoionization 
lifetime and must be continuously replenished (Lecavelier 
des Etangs et al. 1997a). Most interestingly, IUE spectra 
of 51 Oph show variable infalling gas absorption features 
arising from species with a wide range of ionization states, 
including C IV and Si IV, which must be produced by col- 



lisional ionization (Grady & Silvis 1993). Such variable 
features are an indication of the presence of infalling plan- 
ctesimals in the 51 Oph system (for a review of this phe- 
nomenon in the well-studied (3 Pictoris system, see Vidal- 
Madjar et al. (1998)). All this evidence together strongly 
suggests that 51 Oph is similar to j3 Pictoris, and repre- 
sents a young planetary system with a debris disk of gas 
and dust. 

In this paper, we present the first Far Ultraviolet Spec- 
troscopic Explorer (FUSE) spectra of 51 Oph. We observe 
several narrow CS absorption lines arising from the un- 
usual metastable species N I ( 2 D), N I ( 2 P), and S II ( 2 D), 
and calculate the rcdshifts and equivalent widths of these 
lines. We also set low upper limits on the column densities 
of H 2 and CO toward the star, which has interesting conse- 
quences in light of the recent detection of large quantities 
of warm molecular gasses in ISO spectra of 51 Oph (van 
den Ancker et al. 2001). The high S/N FUSE spectra were 
taken 6 days apart, permitting us to detect time- variable 
redshifted absorption from four species, Fe III, N II, N I, 
and metastable O I ( 1 D). We examine the velocity struc- 
ture of the resolved variable features, and carefully deter- 
mine column densities of the infalling gasses. In §7, we 
investigate the unusual composition of the infalling gas, 
and demonstrate that it is most likely produced from in- 
falling planetesimals, as in the case of f3 Pictoris. 

2. OBSERVATIONS AND DATA REDUCTION 



51 Oph was observed with FUSE on 2000 August 29 
and 2000 September 4. The complete spectrum appears in 
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Figure 1. The data were obtained using the low- resolution 
30" x 30" aperture, covering the wavelength range 905 
- 1187 A. The FUSE satellite has four coaligned opti- 
cal channels (LiF 1, SiC 1, LiF 2, SiC 2), forming two 
nearly identical "sides" , each consisting of a LiF grating, 
a SiC grating, and a detector. Each detector is divided 
into two independent segments (a and b), separated by 
a small gap. The entire wavelength range is therefore 
covered with eight partially overlapping spectra that fall 
on different portions of the two detectors (LiF la, LiF lb, 
etc.). A discussion of the on-orbit performance of FUSE 
may be found in Sahnow et al. (2000). The data were cal- 
ibrated with the CALFUSE 1.8.7 pipeline processing soft- 
ware. Detailed information on this software is available at 
http: / / fusc.pha.jhu.edu/analysis / pipeline_reference.html. 
Since none of the observed lines were intrinsically narrow 
enough to use them to determine the true spectral reso- 
lution of our data, we conservatively assume 20 km s _1 
spectral resolution (R = A/ A A = 15000). 

The total exposure time of the first observation was 
9974 s, divided into 23 exposures, and the total expo- 
sure time of the second 9523 s, divided into 20 expo- 
sures. Since the target was not exactly centered within 
the aperture in every exposure, the resulting wavelength 
shifts were corrected by aligning the individual exposures 
before co-addition, using a linear cross-correlation proce- 
dure. Therefore, the image of the star usually fell on a 
slightly different portion of the detector in each exposure, 
greatly reducing the fixed pattern noise in the final co- 
added spectra. In each observation, we obtained the fol- 
lowing signal-to-noise ratios per resolution clement: 69 in 
LiF lb near 1130 A, 20 in SiC 2b near 1080 A, and 36 in 
LiF la near 1050 A. 

The wavelength scales were corrected for the sign error 
in the CALFUSE 1.8.7 heliocentric velocity calculation. 
The uncertainty in the relative wavelength calibration of 
the data is about 5 pixels (~ 9 km s _1 ). In order to deter- 
mine the absolute wavelength calibration of the data, we 
used a rotationally broadened synthetic stellar spectrum, 
overplotted in Figure 1, which included the photospheric 
C I lines seen in the data (Chayer 2001). This synthetic 
spectrum was calculated at the FUSE spectral resolution, 
assuming solar metallicity and the following stellar param- 
eters: T eff = 10000, log g = 4.0, v sin i = 270 km s" 1 . 
The synthetic spectrum does not fit the data well between 
about 1110 A and 1040 A; this may be due to non-LTE 
effects. A paper on improved modeling of the FUV spectra 
of A and late B stars is in preparation. 

In the region of the C I lines, however, the synthetic 
spectrum does fit the data well. The synthetic spectrum 
was shifted to the —21 km s _1 heliocentric velocity of the 
star, then aligned with the co-added LiF lb spectra in the 
region where the photospheric C I lines appear, using the 
cross-correlation procedure mentioned above. This pro- 
cess allowed us to establish the small zero-point offsets of 
the absolute wavelength calibration for the LiF 1 spectra 
with an uncertainty of no more than 14 km s _1 . Since the 
51 Oph spectrum has few photospheric absorption lines 
shortward of 1100 A, the offsets of the absolute wavelength 
calibration for the SiC 2b spectra were determined using 
the narrow IS/CS absorption lines of Ar I. Thus, the un- 
certainty in the absolute calibration for SiC 2b is larger, 



about 27 km s , due to the spread in velocity between 
the unresolved IS and CS components. 

Data taken with FUSE between about 2000 August and 
31 July 2001, when the detector high voltage was raised 
for the second time, may suffer from an effect called de- 
tector x-walk. This effect causes the position of a low 
pulse height photon event to be misplaced in the disper- 
sion direction, and was caused by excess gain sag in local- 
ized places on the detectors where the brightest airglow 
lines fall. The x-walk causes a distortion in the shape of 
an absorption line that conserves photons but not equiva- 
lent width, and varies with time. No good post-aquisition 
remedies for this effect exist at this time. Thus, one must 
be very careful when quantitatively analyzing any absorp- 
tion lines in FUSE spectra taken before the detector high 
voltage was increased, particularly when looking for time- 
variability. To deal with this problem, we carefully exam- 
ined the shape of each line of interest in the spectra from 
two different detector channels, and determined which por- 
tions of the data suffered from the x-walk effect. We ex- 
cluded the few portions of the data that showed significant 
x-walk from our analysis. 

3. OVERVIEW OF FUV SPECTRUM 

The stellar flux from 51 Oph in the FUSE wavelength 
range drops off sharply below 1110 A, as can be seen in 
Figure 1. The most prominent features are the broad C I 
photospheric absorption lines between 1110 A and 1160 A. 
No emission lines, or any other indications of stellar activ- 
ity, are observed in the data. The non-photosphcric ab- 
sorption features in the 51 Oph FUSE spectrum are listed 
in Table 1. Since the absorption features from the lo- 
cal interstellar medium (ISM) are only separated from the 
stellar heliocentric velocity by about 6 km s _1 (Crawford, 
Craig, & Welsh 1997), we are unable to resolve them in 
our spectrum. However, many lines in the spectrum arise 
from excited energy levels, and therefore are not usually 
seen in spectra of the ISM. These excited species indicate 
that there is warm CS gas present near 51 Oph. 

4. METASTABLE SPECIES 

Several of the absorption lines arising from metastable 
states are shown in Figure 2. Only the strongest 
metastable line, arising from O I Q-D), appeared to vary 
significantly between the two days of observation (shown 
in Figure 3c). For the other lines, we averaged the two 
days of data to improve the S/N ratio, and fit Gaussians to 
the lines to determine their redshifts. The strong S II ( 2 D) 
lines were slightly redshifted by 16±14 km s _1 with respect 
to the stellar velocity; the weaker N I ( 2 D) and N I ( 2 P) 
lines appeared to be unshifted. We calculated the equiva- 
lent widths of the lines from the data and converted them 
to column densities, assuming the lines were unsaturated; 
these values appear in Table 2. For the N I ( 2 D 5 / 2 ) lines, 
we had enough lines with different oscillator strengths 
available to determine that the two weakest lines were un- 
saturated, so the column density for this species should be 
accurate. For the N I ( 2 D 3/2 ), N I ( 2 P), and S II ( 2 D) 
lines, this was not the case, so our column densities for 
these species should be regarded as lower limits to the 
true values. In particular, the S II ( 2 D) lines are probably 
considerably saturated. 
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The energies of the metastable states above the ground 
state are 2.38 eV for N I ( 2 D), 3.58 eV for N I ( 2 P), and 
1.84 eV for S II ( 2 D). These energies roughly correspond 
to temperatures of 28000 K, 41000 K, and 21000 K, re- 
spectively. Therefore, these lines must arise solely from 
CS gas; additionally, this gas is quite warm. We calcu- 
lated the radiative decay lifetimes of the metastable states 
by summing the probabilities of the spin-forbidden tran- 
sitions to lower energy levels. These lifetimes appear in 
Table 2. All the lifetimes are quite a bit shorter than the 
six day interval between the two observations; therefore, 
the metastable levels must be continuously repopulated 
by collisional excitation, UV pumping, and/or molecular 
dissociation. Lines arising from the metastable D level 
of C I are seen in spectra of (3 Pictoris (Roberge et al. 
2000) and solar system comets (Tozzi, Feldman, & Festou 
1998). We intend to more fully investigate the significance 
of metastable species in young CS disks in a future paper. 

5. MOLECULAR GASSES 

No molecular hydrogen or carbon monoxide is observed 
toward the star. The 3<r upper limits on the H2 line-of- 
sight column densities in the J = and J = 1 rotational 
levels are N(J = 0) < 3 x 10 13 cnr 2 and N(J = 1) < 
4 x 10 13 cm -2 . The non-detection of absorption in the 
C-X (0,0) band at 1088 A allows us to set an upper limit 
on the line-of-sight column density of 12 CO in the ground 
vibrational level, N( 12 CO) < 7x 10 13 cm" 2 . These low up- 
per limits for the line-of-sight column densities of the two 
most abundant molecular gasses are somewhat confusing 
in light of the recent detection of very large amounts of 
warm molecular gasses in ISO spectra of 51 Oph. Van den 
Anckcr ct al. (2001) detect emission from 12 CO, 13 CO, 
H2O, CO2 and, surprisingly, NO at a gas temperature of 
about 850 K. Our upper limit on N( 12 CO) is highly incon- 
sistent with their column density of 3 x 10 21 cm~ 2 . Even 
at 850 K, such a large column density of 12 CO should 
produce a C-X (0,0) absorption band that is completely 
saturated over several Angstroms. This is clearly not seen 
in our data. Our non-detection of the warm CO seen by 
ISO indicates that very little of it lies along the line of 
sight to the star through the edge-on disk. 

Van den Ancker et al. (2001) propose two possible ex- 
planations for the large amount of molecular gas they de- 
tected, 1) that a large mass ejection from the star into 
the disk has recently occurred, or 2) that a planet-sized 
body orbiting the star has recently been destroyed, pos- 
sibly by a collision with another similar body. Although 
the second scenario seems very exotic, it is more consis- 
tent with our conclusion that very little of the warm CO 
seen by ISO lies along the line of sight. However, this is- 
sue must be thoroughly investigated with more sensitive 
observations searching for absorption from CO and other 
molecular gasses. 

6. VARIABLE FEATURES AND INFALLING GASSES 

Comparison of the spectra from the two days of obser- 
vation shows that they are practically identical. However, 
careful examination of the high S/N spectra revealed a few 
narrow absorption features that varied between the two 
days. They are shown in Figure 3, the N I A1134 triplet in 
Figure 3a, the N II A1085 multiplct in Figure 3b, the line 



arising from metastable Ol ( 1 £>) in Figure 3c, and one of 
the four Fe III lines (1124.87 A) in Figure 3d. No other 
lines in the 51 Oph spectrum showed any significant time 
variation. For comparison, unvarying lines arising from 
C II and ground term O I ( 3 P) arc shown in Figure 4. 
Note that most of the variable lines observed arise from 
excited energy levels, and therefore must represent purely 
CS gas. 

In the August 29 spectrum, all the variable features 
show excess absorption that is redshifted with respect to 
the stellar velocity. The variable infalling gas absorption 
was isolated by dividing the August 29 (day 1) spectrum 
by the quiescent September 4 (day 2) spectrum, removing 
any possible IS, stellar, or non- varying CS absorption; the 
result is shown at the bottom of each panel of Figure 3. 
The spectra from the two days were cross-correlated using 
unvarying lines to align them exactly before normalization, 
but no shifting was necessary to align the LiF 1 spectra. 
All further analysis was performed on these normalized 
spectra. 

The apparent "emission features" in the N II and 
O I ( 1 I?) normalized spectra actually represent excess ab- 
sorption in the day 2 spectrum, in a velocity range that 
is blueshifted relative to the day 1 excess absorption. Al- 
though they appear to be real because they are seen on 
the blue edge of all three N II features, they are not re- 
solved and are only marginally (less than 3er) detected. 
No column densities were calculated for these day 2 excess 
absorption features; data with higher spectral resolution 
would have been necessary to confirm their reality. 

The equivalent widths of the day 1 excess absorption fea- 
tures, calculated from the normalized spectrum, are shown 
in Table 3. They were then converted into column densi- 
ties, assuming that the lines were unsaturated, and us- 
ing the revised absorption oscillator strengths of Morton 
(1999). These column densities represent lower limits to 
the true column densities. More accurate measurement 
of the column densities was performed using the apparent 
optical depth method (see Savage & Sembach (1991) for 
a complete discussion). The apparent optical depth, in 
velocity space, is 



where I (v) is the intensity without absorption, and 
Iobs(v) is the observed intensity. The apparent optical 
depth will be a good representation of the true optical 
depth as long as no unresolved saturated structure is 
present in the absorption line. Unfortunately, for most of 
our lines we cannot prove this directly, as we do not have 
two well separated lines differing in fX arising from the 
same species and energy level. In our case, however, the 
excess absorption extends over at least 40 km s^ 1 , and is 
therefore likely resolved. So, the apparent column density 
in some particular velocity range from V\ to v 2 is 

fV2 fV2 

N a = / N a (v)dv = / r a (v)dv. (2) 

Jvi ne J* JV! 

This method gives accurate column densities for lines that 
are not strongly saturated. A curve of growth analysis 
was not possible, because we did not observe enough lines 
arising from the same species and energy level. Since the 
N II* and N II** lines are blended and may be strongly 
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saturated, we were not able to perform a reliable apparent 
optical depth analysis on them. However, the equivalent 
width analysis at least provides lower limits on the column 
densities of N II* and N II**. 

The velocity ranges of the day 1 excess absorption fea- 
tures and the apparent column densities in these ranges are 
shown in Table 3. First, we can see that many of the lines 
are unsaturated, because the column densities obtained 
from the equivalent widths are equal (within la) to the 
apparent column densities. Second, the column densities 
of the Fe III lines show that the energy levels are statis- 
tically populated, within measurement errors. Therefore, 
the excitation temperature of the infalling Fe III gas must 
be greater than 1340 K, the temperature of the highest en- 
ergy level observed. The presence of O I Q-D) absorption 
implies temperatures of around 23000 K. 

The total column densities of infalling N I and Fe III gas, 
in all fine-structure levels, are A N , = (3.06 ± 0.47) x 10 13 
cm- 2 and A Fe „, = (2.66 ± 0.21) x 10 14 cm" 2 . For N II, 
we find the total column density in all fine-structure levels 
to be at least (2.26 ± 0.17) x 10 14 cm" 2 . The maximum 
possible column density of N II may be calculated by as- 
suming a statistical population of the fine-structure levels 
(1:3:5); the total column density in all levels must be less 
than or equal to (6.7 ±1.0) x 10 14 cm~ 2 . We conclude 
that the total column density of infalling N II gas, in all 
fine-structure levels, is in the range (2.1 — 7.7) x 10 14 cm~ 2 . 

7. DISCUSSION 

Previous workers have suggested that 51 Oph might 
have a companion star, because its radial velocity was 
found to be variable (Buscombe 1963). Our detection 
of time-variable, redshifted CS absorption lines in FUSE 
spectra of 51 Oph confirms the Grady & Silvis (1993) de- 
tection of such lines in IUE spectra. This confirmation 
leads us to believe that any previously measured variation 
in the star's radial velocity is not due to the presence of a 
stellar companion, but rather due to the appearance and 
disappearance of redshifted CS absorption features. 

7.1. Composition of the transient infalling gas 

The ionization energies required to produce Fe III and 
N II are 16.19 eV and 14.53 eV, respectively. Since pho- 
tons with wavelengths less than 766 A and 853 A are re- 
quired to produce these species through photoionization, 
appreciable quantities cannot be produced by the IS UV 
field. Also, as can be seen from the FUSE spectrum of 
51 Oph, the star has virtually no flux below about 1040 
A; therefore, these species cannot be produced by stellar 
photoionization either. They must be produced by colli- 
sional processes involving relatively dense gas. The fact 
that N II is more abundant than N I suggests a moder- 
ately high temperature for the infalling gas. A detailed 
analysis of the ionization balance is beyond the scope of 
this paper, but using the cooling functions of Sutherland 
& Dopita (1993) and assuming collisional ionization equi- 
librium, the ion fraction of N II suggests a temperature 
between 20000 and 34000 K for the transient infalling gas, 
which is consistent with the presence of O I ( 1 .D) in this 
gas. 

The composition of the transient infalling gas appears 
to be rather unusual. No C I, C II, or C III is detected in 



the transient infalling gas. In fact, the only carbon lines 
observed in the FUSE 51 Oph spectrum are photospheric 
C I lines and C II and C II* ground term ( 2 P) lines arising 
from IS/CS gas, which do not vary significantly between 
the two days. The strong metastable C II ( P ) line at 
1010.37 A is not seen, nor are the strong C III lines at 
977.02 and 1175.71 A. Given the estimated gas tempera- 
ture range above, more than 80% of the carbon should be 
present as C II (Sutherland & Dopita 1993). We find a 
conservative 3<r upper limit of 4 x 10 13 cm~ 2 for the col- 
umn densities of both C II and C II* in the infalling gas, 
giving an upper limit on the total carbon column density, 
A c < 1 x 10 14 cm- 2 . 

The total nitrogen column density in the infalling gas, 
N N = An : +A n „, is in the range (2.36— 8.05) x 10 14 cm" 2 . 
Therefore, the ratio of An to Ac in the transient infalling 
gas is at least 2.36. The solar abundance of nitrogen rel- 
ative to carbon, An /Ac, is 0.251 (Grevesse et al. 1990, 
1991). The mean gas phase abundance of carbon atoms 
in the Local ISM is 140 ± 20 per 10 6 H atoms (Cardelli et 
al. 1996) and that of nitrogen is 75 ± 4 per 10 6 H atoms 
(Meyer et al. 1997). These values do not appear to vary 
with direction or with the physical conditions of the IS 
environment. Therefore, the infalling gas in the 51 Oph 
system does not appear to be composed of unprocessed IS 
gas, nor does it have a solar-like composition. 

Given the gas temperature range determined from the 
nitrogen ionization balance, 63% to 64% of the total iron 
should be present as Fe III (Sutherland & Dopita 1993). 
The total column density of iron in the infalling gas, A Fe , 
is therefore roughly (3.83 — 4.56) x 10 14 cm~ 2 . The ra- 
tio of Ap to Ac in the transient infalling gas is at least 
3.83, while the solar abundance of iron relative to carbon 
is 0.0795 (Grevesse & Sauval 1999), indicating that iron is 
even more overabundant relative to carbon than nitrogen 
is. We consolidate the results of our abundance analysis 
below. 

(A N /A C ) 51 0ph > 9 i 

(A N /A C ) S olar " 

(A Fe /A c ) 5 i oph > 4g 2 

(Ap e /Ac) solar " 

(A Fc /A N ) 5 i Qph 
(Ap e /An) solar ^ ' ' 

We were not able to include the oxygen abundance in this 
analysis, because at temperatures of 20000 to 34000 K, 
less than 8% of the oxygen is present as O I. Therefore, we 
require a large extrapolation from our measured O I ^D) 
column density to obtain the total oxygen column den- 
sity. Unfortunately, there are no ground term O II lines 
longward of 911.7 A. 

Natta, Grinin, & Tambovtseva (2000) establish a use- 
ful test to determine if gas infall is related to planetesimal 
evaporation or magnetospheric accretion from a massive 
CS disk, as seen in classical T Tauri stars (Hartman et al. 
1994). If the infalling gas is produced from star-grazing 
planetesimals, it should be depleted in hydrogen (or alter- 
natively, it should be metal- rich). On the other hand, one 
would expect a roughly solar composition if the redshifted 
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absorption is due to magnetospheric accretion. While our 
analysis of the composition of the infalling gas seen with 
FUSE is somewhat incomplete since we have no informa- 
tion on the hydrogen abundance, it does indicate that the 
composition of the transient infalling gas is extremely non- 
solar. Specifically, it appears that the gas must be either 
volatile-depleted or iron-rich, indicating that the gas is 
more likely to have been produced by vaporization of plan- 
etary material. 

7.2. Dynamics of the transient infalling gas 

We calculated the ratio of the force of radiation pres- 
sure from 51 Oph to the force of gravity for each variable 
species observed, using the equations 
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where /j is the absorption oscillator strength, is the 
stellar flux at some distance r from the star, M sta r is the 
mass of 51 Oph, and M at0 m is the mass of the atomic 
species. Note that the ratio of the forces, (3, is inde- 
pendent of the distance r. The stellar fluxes were either 
measured from the FUSE spectrum, or estimated from 
our synthetic stellar spectrum at longer wavelengths. The 
strongest transitions are in the FUSE wavelength range, 
and the values of (3 do not depend sensitively on the syn- 
thetic spectrum fluxes. We used a stellar mass of 3 M Q , 
estimated from the star's position on the H-R diagram and 
assuming that it is on the zero-age main sequence. The cal- 
culated values of (3 are shown in Table 3. We can see that 
for most of the variable species (3 is greater than 1, and 
for N I and O I ( 1 -D) it is quite large. One would therefore 
generally expect any optically-thin absorption from these 
species to be blueshifted, unless the gas was injected to- 
ward the star at a relatively high radial velocity, as would 
happen if it was produced from infalling planetesimals. 

8. CONCLUDING REMARKS 



A comparison of the characteristics of 51 Oph to those of 
(3 Pictoris shows many similarities, strengthening the sug- 
gested evolutionary connection between the two systems. 
Both CS disks contain highly excited, metastable species 
that have short radiative decay lifetimes. Both systems 
also show the presence of variable, collisionally ionized CS 
material. The variable absorption features in both systems 
are most frequently redshifted with respect to the stellar 
velocity, despite the effect of radiation pressure. There are 
some differences between the two systems, however. If the 
variable redshifted events in the 51 Oph system are caused 
by infalling planetesimals, this system has a lower infall 
rate than (3 Pic, since many observations of 51 Oph have 
not shown any variable infalling gas features (Lecavelier 
des Etangs et al. 1997b; Crawford, Craig, & Welsh 1997). 

The composition of the transient infalling gas in 51 Oph 
is highly non-solar, and is more likely produced by the 
destruction of a star-grazing planetesimal than by mag- 
netospheric accretion of unprocessed gas. However, this 
conclusion needs to be tested by observation of additional 
infalling gas events in lines arising from many species, in 
order to perform a more complete analysis of the infalling 
gas composition. And detailed modeling is necessary to 
fully understand the significance of the metastable atomic 
species in both (3 Pictoris and 51 Oph. Identification of 
51 Oph as a planetary, rather than protoplanetary, sys- 
tem will at last firmly add another object to the class of 
(3 Pictoris-type systems. Since 51 Oph is much younger 
than [3 Pictoris, which is about 2 x 10 7 yrs old (Barrado 
y Navascues et al. 1999), this will provide a tighter con- 
straint on the time scale for planetary formation, and in- 
crease our knowledge about the evolution of CS disks. 
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Table 1 

is and cs absorption lines in the fuse spectrum of 51 oph 



Species 


Wavelength 


Energy Level 


Comments 


(A) 


(cm" 1 ) 




CII 


1036.34 


0.00 


part IS 


CII* 


1037.02 


63.42 


purely CS 


OI 


1039.23 


0.00 


part IS 


OI* 


1040.94 


158.27 


purely CS 


O I** 


1041.69 


226.98 




Si II 


989.87, 1020.70 


0.00 


part IS 


Si II* 


992.69, 1023.70 


287.24 


purely CS 


Arl 


1048.22, 1066.66 


0.00 


part IS 


Fell 


1142.37, 1055.26 


0.00 


part IS 


Fell* 


1146.95, 1060.44 


384.79 


purely CS 


Fe II** 


1150.29, 1063.02 


667.68 


5) 


Fe II*** 


1153.27, 1065.84 


862.61 


5> 


Fe II**** 


1154.40, 1076.85 


977.05 





Metastable Species 



N I 
N I 
N I 
SII 
S II 



( 2 D 5/2 ) 
( 2 D 3/2 ) 

CD,/ 

( 2 D : 



2) 

5/2) 



1163.88, 1176.51 

1168.54, 1177.70 
1143.646, 1143.651 

1019.53 

1014.44 



19224.47 
19233.18 

28839 
14852.87 
14884.67 



Time-Variable Species 



purely CS 
purely CS 



N I 


1134.17, 1134.41, 1134.98 


0.00 


part IS 


N II 


1083.99 


0.00 




N II* 


1084.57, 1084.58 


48.67 


purely CS 


N II** 


1085.53, 1085.55, 1085.71 


130.80 




OI CD) 


1152.15 


15867.86 


purely CS, metastable 


FeUI 


1122.52 


0.00 




Fe III* 


1124.87 


435.80 


purely CS 


Fe III** 


1126.72 


738.55 




Fe III*** 


1129.19 


932.06 
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Table 2 
Metastable CS SPECIES 



Species 


Wavelength 

(A) 


Energy level 
(cm- 1 ) 


Lifetime 


(mA) 


N Wx h 
(cm- 2 ) 


Comments 


N I ( 2 D 5/2 ) 
N I ( 2 D 3/2 ) 
N I ( 2 P) 
S II ( 2 -D 3/2 ) 
S II ( 2 -Db/ 2 ) 


1168.42 
1177.70 
1143.646, 1143.651 
1019.53 
1014.44 


19224.47 
19233.18 

28839 
14852.87 
14884.67 


48.2 hrs 

12.3 hrs 
13 sec 
39 min 
36 min 


2.5 ± 1.1 
8.4 ± 1.3 
7.98 ±0.92 
54.0 ± 7.6 
86.7 ±5.8 


(1.06 ±0.45) x 10 14 
(5.07 ±0.78) x 10 13 
(5.77 ±0.66) x 10 12 
(1.15 ±0.16) x 10 15 
(1.68 ±0.11) x 10 15 


unshifted 
unshifted 
unshifted, blended a 

shifted 13 

shifted 13 



a For the blended lines, we used an effective oscillator strength equal to the sum of the two individual line oscillator 
strengths. 

b The metastable S II lines are slightly redshifted by 16 ± 14 km s _1 with respect to the stellar velocity. 



Table 3 
Variable infalling gasses 



Species 


Wavelength 


Energy level 




Nw x b 


Velocity range c 


/3 d 


N a c 


Comments 




(A) 


(cm- 1 ) 


(mA) 


(10 13 cm" 2 ) 


(km s- 1 ) 




(10 13 cm" 2 ) 




N I 


1134.41 


0.00 


9.1 ±2.8 


2.68 ±0.83 


0-40 


7.4 


3.07 ±0.76 


unsaturated 


N I 


1134.98 


15 


15.7 ±3.2 


3.16 ±0.65 


15 


15 


3.06 ±0.47 




N II 


1083.99 


0.00 


58 ± 13 


4.8 ± 1.1 


14 - 114 


1.6 


7.4 ± 1.2 


saturated 


N II* 


1084.57, 1084.58 


48.67 


85 ± 11 


7.09 ±0.91 




1.6 




blended 1 


N II** 


1085.53, 1085.55, 1085.71 


130.80 


95 ± 11 


8.02 ±0.92 




1.6 




blended 1 


OI CD) s 


1152.15 


15867.86 


22.7 ±3.5 


1.84 ±0.28 


6-71 


4.7 


1.85 ±0.30 


unsaturated 


FeUI 


1122.52 


0.00 


48.5 ±4.4 


7.99 ± 0.73 


0-84 


0.25 


9.42 ±0.85 


saturated 


Fe III* 


1124.87 


435.80 


26.4 ±2.7 


6.58 ±0.67 


15 


1.3 


6.26 ±0.68 


unsaturated 


Fe III** 


1126.72 


738.55 


11.4 ±3.1 


5.4 ± 1.5 


1) 


0.68 


6.6 ± 1.5 




Fe III*** 


1129.19 


932.06 


16.9 ±4.3 


3.75 ±0.96 


55 


0.46 


4.32 ± 0.97 





a Equivalent width. 

b Column density calculated from W\, assuming the line is unsaturated. These values are lower limits to the true column densities. 
c Velocities given with respect to the stellar velocity. 
d Ratio of radiation pressure to gravity for this species. 

° Apparent column density, obtained by integrating the apparent optical depth over the velocity range. 

f For the blended lines, we used an effective oscillator strength equal to the sum of the individual line oscillator strengths. 

s Metastable species. 
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Fig. 1. — Overview of the entire 51 Oph FUSE spectrum. The data from the two days of observation have been averaged together and 
rebinned by a factor of 32 for this plot. The synthetic stellar spectrum described in the text is shown in red. The prominent C I photospheric 
absorption lines between 1110 A and 1160 A are indicated. The Lyman-/3 airglow line is marked with ffi. 
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Fig. 2. — Metastable species in the 51 Oph CS gas. The ± la measurement error bars are overplotted in red. In the top two panels, data 
taken 2000 August 29 are plotted in black and those taken 2000 September 4 in green. None of the lines showed any significant variation 
between the two days of observation. In the bottom two panels, the data from the two days of observation were averaged together to improve 
the S/N ratio, (a) Two of the N I ( 2 D) lines in the LiF lb segment spectra, rebinned by a factor of 2. (b) The blended N I ( 2 P) A1143.65 
doublet in the LiF lb segment spectra, rebinned by a factor of 2. The FUSE spectral resolution is insufficient to separate the two lines of the 
doublet. (c) The S II ( 2 -D 5/2 ) A1014.44 line in the LiF la segment spectrum, rebinned by a factor of 4. (d) The S II ( 2 D a/2 ) A1019.53 
line in the LiF la segment spectrum, rebinned by a factor of 4. 
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Fig. 3. — Time-variable absorption lines in the 51 Oph FUSE spectra. Data taken 2000 August 29 (black) and 2000 September 4 (green) 
are shown at the top of each panel. The August 29 spectrum divided by the September 4 spectrum is shown at the bottom of each panel. 
The ± la measurement error bars are overplotted in red. Velocities are given in the rest frame of the star, (a) The N I A1135 triplet in the 
LiF lb segment spectra, rebinned by a factor of 2 for this plot. (b) The N II A1085 multiplet in the SiC 2b segment spectra, rebinned by a 
factor of 8. (c) The O I ( x O) A1152.15 line in the LiF lb segment spectra, rebinned by a factor of 2. (d) The Fe III* A1124.87 line in the 
LiF lb segment spectra, rebinned by a factor of 2. 
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Fig. 4. — Unvarying absorption lines in the 51 Oph LiF la spectra, showing data taken 2000 August 29 (black) and 2000 September 4 
(green). No variable infalling gas is seen in these lines arising from carbon and oxygen. The data were rebinned by a factor of 4 for this plot. 
The ± 1<t measurement error bars are overplotted in red. 



